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We investigate the influence of the central two-body nucleon-nucleon correlations on several quan- 
tities observed in relativistic heavy-ion collisions. It is demonstrated with explicit Monte Carlo 
simulations, that the basic correlation measures observed in relativistic heavy-ion collisions, such as 
the fluctuations of participant eccentricity, initial size fluctuations, or the fluctuations of the number 
of sources producing particles, are all sensitive to the inclusion of the two-body correlations. The 
effect is at the level of about 10-20%. Moreover, the realistic (Gaussian) correlation function gives 
indistinguishable results from the hard-core repulsion, with the expulsion distance set to 0.9 fm. 
Thus, we verify that for investigations of the considered correlation measures, it is sufficient to use 
the Monte Carlo generators accounting for the hard-core repulsion. 
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I. INTRODUCTION 

The atomic nucleus is closer to a self-bound saturated 
liquid than to a Fermi gas of non-interacting particles, 
as is for simplicity frequently assumed in studies of rel- 
ativistic heavy-ion collisions. Thus the inclusion of cor- 
relations in the initial configuration of nucleons in the 
colliding nuclei is a priori very important. Recently Alvi- 
oli, Drescher, and Strikman generated distributions 
of nucleons in nuclei vi^hich account for the central two- 
body nucleon-nucleon (NN) correlations. The procedure, 
based on the Metropolis search for configurations satis- 
fying constraints imposed by the NN correlations, repro- 
duces the one-body Woods-Saxon distributions, as well 
as central NN correlations, taken in the Gaussian form. 
This calculation is a very important step in the investi- 
gations using the Glauber approach dJil to relativistic 
heavy- ion collisions, as it is well known |5|,|6| that correla- 
tions induce event-by-event fluctuations of the measured 
quantities. 

The Glauber Monte Carlo codes 

(for a discus- 
sion of physics issues see Ref. [§] and the review (TT| ) 
which model the early phase of the collision, have not 
been incorporating, for practical reasons, realistic NN 
correlations. Instead, the hard-core expulsion, easy to 
implement, is used. In that method, centers of nucleons, 
whose positions are randomly generated according to the 
Woods-Saxon one-body distribution, are not allowed to 
be placed closer to one-another than the expulsion dis- 
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tance d ^ 1 fm, which simulates the hard-core NN re- 
pulsion. It is not a priori clear that the results obtained 
with the realistic (Gaussian) and the hard-core correla- 
tions should be the same for various correlation measures 
used in the heavy-ion studies. Moreover, it is not obvi- 
ous what precise value of d should be taken to make the 
simulations most realistic. 

The purpose of this paper is to investigate, with the 
help of explicit Glauber Monte-Carlo simulations by 
GLISSANDO the role of the central two-body NN cor- 
relations for several popular observables in relativistic 
heavy-ion collisions. In particular, we look at the fol- 
lowing fluctuation measures: the participant eccentric- 
ity fluctuations related to the fluctuations of the elliptic 
flow p"3 - [2^ . the multiplicity fluctuations as analyzed in 
the set up of the CERN NA49 experiment |25| . and the 
recently investigated initial size fluctuations [2a | , w hich 
influence the transverse-momentum fluctuations |27h42 | . 
We find that all these measures are sensitive to the in- 
clusion of the two-body correlations at a level of about 
10-20%. However, the realistic (Gaussian) correlation 
function gives virtually indistinguishable results from the 
calculations with the hard-core repulsion, with the expul- 
sion distance tuned to c? = 0.9 fm. Thus, we will argue 
that for all practical terms of modeling the Glauber ini- 
tial phase of the collision, it is sufficient to use the Monte 
Carlo generators with the hard-core repulsion. 

Certainly, the method of Ref. ^ is more general, as 
it allows to include correlations from attractive forces, 
as well as introduce the isospin dependence. These were 
recently considered in Ref. |4j] , and when these distribu- 
tions are published, they can be implemented in Glauber 
generators and tested in a similar way as in the present 
work. 



2 




U 



1.0 

0.8 - 

0.6 - 

0.4 - 

0.2 - 
0.0 - 



'i 



\ 



0.25 
7 0.20 
^ 0.15 
"q. 0.10 
^ 0.05 

0.00 

r [fm] 

FIG. 1: (color online). Radial one-body density of centers 
of nucleons in ^"^Pb, A-Kr^ p^^Ur) , obtained numerically from 
the distributions of Ref. [J, (dots). The line shows the 
Woods-Saxon fit with the optimum parameters ((2]). 
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FIG. 2: (color online). Central two-body NN radial correla- 
tion density for ^"^Pb, obtained from Eq. ([3} (points), and 
the Gaussian fit of Eq. (|5I6I) (line). 



II. NUCLEAR CORRELATIONS 

The method of Ref. [l| imposes a given form of one- and 
two-body nucleon distributions. The one-body density is 
parametrized with the standard Woods-Saxon form 



p(i)(r) = 



A 



(1) 



Our fit to the distributions for ^°®Pb from [5] yields the 
optimum parameters 



R = 6.59(1) fm, a = 0.549(2) fm. 



(2) 



where the uncertainties follow from the regression analy- 
sis on the available sample [2] of 10^ configurations. The 
result of our numerical simulation is displayed in Fig. [TJ 
The radial two-body correlation function C(r) is de- 
fined as [l| 



Pu in 



(3) 



f2l (2) 

where Pq (r) and Pjj \r) are the correlated and uncorre- 
lated radial two-body densities. 



Pf\r) 



Here p\ (ri, r^), i — C,U, denotes the appropriate two- 
nucleon density, r is the relative coordinate, r — |r|, and 
il corresponds to the two angles associated with r, over 
which the density is integrated. The correlated density 
is read off from the distributions [2I with the help of 
GLISSANDO by histogramming the relative distances be- 
tween the centers of nucleons in the same nucleus, while 
the uncorrelated density is found by taking the pairs of 
nucleons from different nuclei (this corresponds to the 
well-known mixing technique, which gets rid of correla- 
tions). The result of our procedure is shown in Fig. [21 



d'n I d^RpP{R + r/2,R-r/2). (4) 



We recover the Gaussian central NN correlation, imple- 
mented in the procedure of Ref. pj , 



C{r) 



e 2b-i 



with 



5 = 0.561(1). 



(5) 



(6) 



The uncertainty comes from the finite sample of 10"" con- 
figurations from 

Thus indeed the distributions of [T, properly im- 
plement the one-body density and the Gaussian central 
two-body correlations. The purpose of the above study 
was to read off the one-body parameters which in 
the following sections will be input in the generation of 
the uncorrelated distributions by the Glauber simulations 
with GLISSANDO Results from the uncorrelated dis- 
tributions will be compared to the correlated case, where 
the correlated distributions of Ref. fT] will be fed directly 
into our simulations. 



III. GLAUBER MODELS 

The prototype Glauber model used in the h eavy -ion 
phenomenology is the wounded-nucleon model |44l |. A 
wounded nucleon has collided inelastically at least once 
in the collision process. Variants of the approach l45l - 
l47j admix a certain fraction of binary collisions to the 
wounded nucleons, which leads to a better overall de- 
scription of multiplicities of the produced particles. In 
this mixed model, investigated in this work, the number 
of the produced particles is proportional to the number 
of sources 



(1 - a)iV„,/2 + a7Vbi, 



(7) 



where is the number of the wounded nucleons, and 
A'bin the number of binary NN collisions. The fits to 
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particle multiplicities of Ref. [47[ give a = 0.145 for 
collisions at ^/snn = 200 GeV, and a = 0.12 for 
y'SNN = 19.6 GeV. Extrapolation to the LHC energies 
yields a ~ 0.2. 

More sophisticated approaches psi - lsoj discriminate 
between the nucleons which have collided only once 
(corona) and more than once (core), which leads to an ap- 
pealing physical picture. Also, the wounded-quark model 
[sH - ISq yields a successful phenomenology. All in all, the 
Glauber picture of the initial stage of the relativistic 
heavy-ion collision is a key element of many phenomeno- 
logical analyses of the particle production mechanism. 

In this paper we apply the mixed model for the ^^^Pb- 
208pj-, collisions, with a = 0.12, corresponding to the 
highest SPS energy. We term the locations of cen- 
ters of the wounded nucleons or the binary collisions 
as "sources", with the weight of the wounded nucleon 
Wi — (1 — a)/2, and the weight of the binary collision 
Wi = a. A source emits particles, according to a super- 
posed distribution 0. 

While for the one-body measures, such as the particle 
multiplicities or spectra, only the one-body distributions 
matter and correlations are irrelevant, the fluctuations 
measures are expected to be sensitive to the NN corre- 
lations in the nucleon distributions. These are examined 
in detail in the next section. 



IV. RESULTS OF SIMULATIONS 

In this section we compare the results of the Glauber 
calculation initialized with the distributions of Ref. [l|, ^ 
(solid lines in the figures), with uncorrelated distributions 
(dashed lines), and with the distributions accounting for 
the hard-core repulsion with the expulsion radius d = 
0.9 fm (dotted lines). The simulations are performed 
with GLISSANDO ^. 

We note that in the case with no correlations we sim- 
ply use the Woods-Saxon parameters ([2]), while in the 
case with the hard-core repulsion we need to start with 
a somewhat more compact distribution, as the expulsion 
leads to swelling, as explained in Ref. Q. We find that 
starting the Monte Carlo generation with R = 6.44 fm 
and a = 0.549 fm, leads, with d — 0.9 fm, to the one- 
body distribution with parameter values This con- 
struction, with shrunk "bare" one-body distributions, is 
important, as that way all calculations presented in the 
figures correspond to identical one-body distribution, and 
the differences in results are caused entirely by the two- 
body correlations. 



A. Eccentricity 

We start with a measure sensitive to the fluctuations, 
the so-called participant eccentricity. This measure ap- 
pears in the studies of the event-by-event fluctuations 



of the initial shape, in particular of its elliptic compo- 
nent [T^ - [2^ . The eff'ect is important, as the fluctuations 
lead to enhanced eccentricity of the initial system, and 
as a result of the subsequent hydrodynamic evolution, 
to enhanced elliptic flow. The participant eccentricity is 
defined in each event as 



fo-z _ 0-2)2 _^4o-2 



(8) 



where cr^ ^nd cr^ are the variances of the two transverse 

X y 

coordinates, and a^y is the covariance. Specifically, in 
each event 



{x) = ^ WiXi, = ^ w^{xi - {x)Y 



(9) 



and similarly for the y variable and the covariance. The 
index i runs over all generated sources, and Wi are the 
weights. The quantity e* has the interpretation of the ec- 
centricity evaluated event- by- event in a variable reference 
frame jl9i | , rotated in such a way that the eccentricity in 
a given event is maximized. 

In the top panel of Fig. [3] we show the dependence 
of the event-by-event average, (e*), on the number of 
wounded nucleons (determining the centrality of the 
event). We note that the three calculations are virtu- 
ally indistinguishable, except for a tiny difference for the 
most central collisions, where the uncorrelated case is a 
few percent higher. The same conclusions were reached 
in the analogous study of eccentricity in Ref. [57| . 

The bottom panel of Fig. [3] shows the scaled standard 
deviation, Ae*/(e*), obtained from our event- by- event 
analysis. We note a significant difference between the 
uncorrelated case, which has up to 10% larger fluctu- 
ations at intermediate centralities, and the cases with 
correlations. However, the calculations with the realistic 
NN correlations and the hard-core correlations give an 
indistinguishable result, with the two curves overlapping 
within the statistical noise. 

The short horizontal line at the most central events 
corresponds to the theoretical value ^4/7r — 1 of 
Ref. [1^, following from the central limit theorem. 



B. Multiplicity fluctuations 

Next, we consider a quantity relevant for the multi- 
plicity fluctuations as measured in the NA49 experimen- 
tal setup .25.] . where the number of participants in the 
projectile is determined via the VETO calorimeter. Sig- 
niflcant fluctuations of the number of sources may follow 
in this case from the fact that even at a flxed number 
of the wounded nucleons in the projectile, the number of 
wounded nucleons in the target fluctuates due to the sta- 
tistical nature of the Glauber approach. The fluctuations 
of multiplicity in nucleus-nucleus collisions were also in- 
vestigated experimentally in [58l - l6l| . We recall [H, US] 
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FIG. 3: (color online). Top: the average participant eccentric- 
ity, (e*), vs the number of wounded nucleons, N^^, obtained 
with the three investigated nucleon distributions described in 
the text. Bottom: The scaled standard deviation Ae*/(e*) 
obtained from an event-by-event study. The short horizontal 
line at the most central events corresponds to the theoretical 
value \f^jTF~- \ of Ref. [l^ following from the central limit 
theorem. 
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FIG. 4: (color online). The scaled variance of the number of 
sources, a;, plotted as a function of the wounded nucleons in 
the projectile, iV^^°-'. 
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FIG. 5: (color online). The scaled standard deviation of the 
size variable r of Eq. (|lip . plotted as a function of the total 
number of wounded nucleons, A^^,. 



that the simple superposition models with the effect of 
fluctuations of the target wounded nucleons are not able 
to explain the data of Ref. [1^. Nevertheless, for the 
present purpose of analyzing the importance of the NN 
correlations, the effect serves its purpose. 

In Fig.|4]we show the scaled variance of the total num- 
ber of sources defined in Eq. (O , 



var(iV,) 



(10) 



plotted as a function of the wounded nucleons in the 
projectile, N^^'^ . We note a significant, about 20%, 
reduction of w when the two-body NN correlations are 
included. However, again there is no noticeable differ- 
ence between the realistic (Gaussian) correlations and 
the hard-core expulsion, as the two lower curves in the 
figure overlap. 



C. Size fluctuations 

Finally, we look at the event-by-event size fluctuations, 
namely the fluctuations of the variable 



(11) 



It was recently shown in Ref. [26'] that the initial size 
fluctuations are carried over via hydrodynamics and sta- 
tistical hadronization into the event-by-event transverse- 
momentum fluctuations [27l - l4^ . where they lead to a 
natural description of the RHIC data for the measure 
o'dyn(PT)- In Fig.[S]we show the scaled standard deviation 
of r. Once again, the presence of the NN correlations re- 
duces somewhat the fluctuations, while the realistic and 
hard-core correlations with d — 0.9 fm give virtually the 
same result. 
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V. CONCLUSIONS 

We have checked by carrying out exphcit Glauber 
Monte Carlo simulations with GLISSANDO Q, that the 
inclusion of the central NN correlations influences the 
fluctuation measures in relativistic heavy-ion collisions 
at a level of, say, fO-20%. Comparison of the realis- 
tic (Gaussian) correlations implemented in Ref. [l| and 
the hard-core correlations, typically used in the Glauber 
Monte Carlo codes, shows that they lead to the same re- 
sults when the hard-core expulsion distance between the 
centers of nuclcons is tuned to 

d = 0.9 fm. (12) 

Thus the main massage for the practitioners of the 
Glauber Monte Carlo models is that, at least for the in- 
vestigated observables, the implementation of the hard- 



core repulsion with d given by Eq. p2|) . straightforward 
to implement in Monte Carlo generators, leads to realistic 
predictions. We note that the dependence of the results 
on the value of d is sensitive, as the excluded volume 
scales as d^. 



Certainly, the general method of Ref. [l| allows one to 
implement channel-dependent NN correlations, as well as 
the nuclear attraction, relevant at intermediate distances. 
The role of these effects for the fluctuation measures in 
relativistic heavy-ion collisions can be investigated in a 
similar manner as in this work. In essence, every effect 
which increases the "regularity" of the initial nucleon dis- 
tributions of the colliding nuclei, such as the considered 
central NN correlations, will have the tendency of de- 
creasing the event-by-event fluctuations in nuclear colli- 
sions generated by the Glauber models. 
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